ADVANCED VISION INTERVENTION ALOGRITHM 

1 Cross Reference to Related Application : 

2 This application is a continuation of Provisional Patent Application Serial No. 60/250,296, 

3 filed November 30, 2000. 

4 Field of the Invention: 

5 This invention relates in general to the field of optically corrected aberrant human vision and 

6 more particularly, relates to the construct of an improved optical model of the human eye and 

7 procedures to achieve a more precise optical correction. 
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p Background of the Invention : 

Li 

j}. The human eye often contains aberrations that may be optically corrected for normal 

ip activities. If the primary aberration causes a distant object to be focused in front of the retina, for 

L = 
r~ 

\'U example, the eye is said to be myopic or near-sighted. If the primary aberration causes the distant 
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llf object to be focused behind the retina, the eye is said to be hyperopic or far-sighted. Often the 

13** aberrations are such that the light rays which enter the cornea in one location from a distant object 

14 focus at one point inside the eye while light rays which enter the cornea in another location from the 

15 same distant object focus at a different point inside the eye. A typical example of this is called 

1 6 astigmatism in which light rays which enter the eye in one meridian focus at one point and light rays 

1 7 which enter in the meridian perpendicular focus at another location. It is easy to visualize that other 

18 more general aberrations can and do occur in the human eye. 

19 Interventions in common use today to correct such human eye aberrations are spectacles, 
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contact tenses, corneal refractive surgery, such as LASIK or corneal implants, and phakic intraocular 
lenses. The correct specification of spherocylindrical spectacles and contact lenses is well 
established. The newer surgery based techniques are not as predictable as may be desired, but still 
are in wide spread use. 

A current development is instruments intended to measure the total aberrations of the eye 
rather than simply the traditional spectacle correction parameters of sphere, cylinder, and axis. These 
instruments provide the theoretical information required to improve on vision correction beyond the 
lower-order aberrations of sphere, cylinder, and axis. The application of this information in a clinical 
setting has not been entirely successful because the optical measurements of the eye and the 
intervention methods are not exact. 

Another area of interest is the relatively mature process of calculating the intraocular lens 
(IOL) power for post cataract patients. The calculation of implantable lens powers and prediction 
of post-operative outcomes has been an ongoing research interest since the first implanted lens. The 
calculations employed by modern strategies may effectively reduce the prediction error bias for the 
majority of normal cases but even with these adaptive third-generation IOL calculations, there are 
a significant number of real world cases whose outcomes are poorly predicted. A particular category 
that is poorly predicted is cataract surgery following previous refractive surgery. A primary source 
of error in the corneal power measurement with keratometers is that these instruments typically 
measure corneal curvature in the 32 mm zone and can miss the relatively flat central region. Modern 
IOL calculation strategy, as stated above, is based on calculations involving a mix of parameters 
intended to improve the effective lens position predictor (ELP) and calculation. Attempts to improve 
ELP have employed such parameters as axial length, eomeal power, horizontal corneal diameter, 



anterior chamber depth, crystalline lens thickness, preoperative refraction, and age. The basic 
methodology for the modern IOL calculation strategy has three steps: a) compute the ideal theoretical 
power of the implanted lens (this lens power will not necessarily be available) using measured and 
adapted optical model parameters, b) for the surgeon selected available power near the ideal power, 
predict the post operative refractive outcome, and c) given the actual outcome of the procedure, 
update optical model parameters to decrease the prediction error of future surgeries. As stated above, 
such modern strategies, while improving the prediction error bias for the majority of normal cases, 
do not improve predicability for many other cases, such as the case of cataract surgery following 
previous refractive surgery. One modern attempt at improvement is custom LASIK using corneal 
topography and/or wavefront aberration data. However, large amounts of positive spherical 
aberration can be induced as a result of LASIK procedures and for large pupil diameters, the retinal 
image quality may be degraded. 

While past and current developments have improved the predictability of optical correction 
of aberrations in the human eye, as may be understood by those skilled in the art, the presently 
known techniques or methodologies have not produced a measure of predictability for the many 
variations found in the parameters describing the human eye. 



Summary of the Invention : 

The inventive principles of the advanced vision intervention algorithm (A VIA), as disclosed 
in the description of the preferred embodiment, may be used to improve the predictability of most 
known methodolgies for optical vision correction and may be applied to any similar methodology 

\ 

A? which may be developed. The advanced vision intervention algorithm (A VIA) provides a means to 



1 optimize the predictability of almost any current or anticipated customized human vision 

2 intervention method. In the following summary and Description of a Preferred Embodiment, AVIA 

3 is described in general terms in a series of steps. For the purpose of explanation, the AVIA method 

4 is shown in nine steps but as would be known by those skilled in the art, the number of steps used 

5 to describe the invention may be varied without departing from the principles of the disclosed 

6 invention. 

7 For the example shown in the preferred embodiment, a summary of these steps is: 

8 1. Input data: The input data is evaluated for obvious errors. Individual classes are used 

9 for each data type so that new data types and their evaluation methods are easily added to the 
10| framework. This includes inputting a set of categorical data for an eye and visual correction method 

and inputting a set of continuous data for an eye and visual correction method. The input data 

Ij2_ expresses the patient's visual correction need in terms of wavefront aberration. If the wavefront 

it" aberration is not explicitly measured using a system intended for this purpose, it may be estimated 

Li 

lpte using the information provided such as spectacle 

ill 2. Build Optical Model: The second step is building an optical model of the eye and 

W" visual correction method that has been adapted based on prior visual correction outcomes. This 

1 7 processing step employs a ray transfer element (RTE) that is described in detail below. The RTE is 

1 8 an important component of AVIA. 

19 3. Generate Initial intervention recommendation. In this step an initial intervention 

20 recommendation is computed based upon the optical model. This recommendation will either be 

21 categorical (e.g., selection from a series of available implantable lenses); a continuous scalar such 

22 as a theoretical power value; or a continuous vector quantity such as a surface description, or a 

4 



combination of these such as a lens series with a discrete set of base curves and a customizable front 
surface. \ 

4. Outcome Prediction: In this step, the outcome of the intervention is predicted in terms 
of the wavefront aberration^ The predicted wavefront aberration is displayed in a number of ways 
any one of which can be selected by the surgeon for viewing. Examples of these displays include 
equivalent spectacle lorrection, wavefront aberration variance, wavefront aberration contour map, 
point spread function metrics and displays, modulation transfer function metrics and displays, and 
simulated retinal images. 

5. Predicted Outcome Evaluation: In this step the surgeon uses the various displays of 
the predict ed outcome waVefront aberration to decide if the outcome would be acceptable for the 
current case under consideration. 

6. Intervention Design Iteration: In this step, the surgeon may continue iterating a 
modification of the Intervention plan and re-evaluates the predicted outcome until the predicted 
outcome is judged as acceptable. 

7. Perform Intervention : In this step the surgeon performs the intervention. 

8. Evaluation of Outcome: In this step the outcome of the intervention is assessed using 
data such as wavefront aberration exams, uncorrected and best corrected visual acuity, contrast 
sensitivity test, post-intervention spectaeie correction, etc. 

9. Update Historical Database\ In this step the database of historical visual correction 
cases is updated to include the current case\ 



1 Brief Description Of The Drawings : 

2 Figure 1 shows in a block element flow diagram, the method of a preferred embodiment, 

3 according to the disclosed inventive principles; 

4 Figure 2 shows in a simplified optical system for the human eye, how step 2, as described 

5 in the Summary and Description of a Preferred Embodiment and shown as numeral 2, in Figure 1, 

6 for the modeling and use of a RTE; 

7 Figure 3 shows an intersection point of an incident ray 1, and the RTE referenced to the 

8 meridional and transverse planes; and 

9 Figure 4 shows the basis vectors for the meridional coordinate system of Figure 3. 
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Description of a Preferred Embodiment : 

N 

lS? The invention, according to its inventive principles, is disclosed with reference to the 

1=3 preferred embodiment, as described below, and with reference to the accompanying drawings. As 

< 

lk* shown in step 1 in Figure 1 , the first step uses the clinical and categorical data of the eye to construct 

1§;J an optical model, as shown in Figure 2, and in step 2 of Figure 1 . 

\ r 6 Optical Model 

17 The foundation of the optical calculations and predicted outcome of a given intervention 

1 8 strategy is the optical model. The optical model is created using all available data for a given eye. 

1 9 This data could be as simple as the two basic quantities of keratometric data and spectacle correction 

20 or as complex as specifying multiple surfaces, axial locations, and wavefront aberration data, or 

21 additional input data, without departing from the disclosed inventive principles. To explain the 

22 procedure, a simple model consisting of a corneal surface and wavefront aberrations only, is used. 
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In Figure 2 is shown the cornea and wavefront W positioned along an optical axis. In this Figure 2, 
the wavefront is located at the entrance pupil of the eye. As shown in Figure 2, if a ray R2 is traced 
from the fovea, it would appear to exit the eye in the direction of RO. This exiting ray is 
perpendicular to W as shown by the dotted line intersection with W. 

To constSruct a model consistent with this information an average axial length of 24 mm may 
be assumed and a ray transfer element (RTE) introduced that maps an incoming ray R2 from the 
fovea to the intermediate ray Rl so that the exit ray RO has the desired orthogonality relationship 
with W, as shown in Figure 2. 

Constructiomof the RTE is accomplished using the following steps: 

1 . For a regular sampling of points on W do the following steps: 

2. Construct theVay RO perpendicular to W at the sample point on W; 

3. Trace - RO through the system until it intersects the RTE plane; 

4. At the RTE intersection point, compute the mapping for the incident and transmitted rays; 
and \ 

5. If there are more samples on W go back to step 2. 

The optical ray tracing of refractive and reflective elements is well known to those skilled 
in the art, for example as shown in JgfJd cm Optical Engin e ering by W.J. Smith an er-pnblished by 
McGraw"Hilh 19 9 0 , and for that reasori is not described in detail. 
Ray Transfer Element \ 

The ray transfer element (RTE) is an important aspect of AVIA. It provides both a means to 
generate an optical model that is consistent witttvexam data and a means to adapt the modeling to 
past surgical procedures so that prediction of postoperative outcomes is improved for future 
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surgeries. The basic operation of the RTE is to efficiently map incident rays to transferred rays 
(either transmitted or reflected rays) and provide a method to smoothly adjust the mapping between 
data points to simulate continuous data. Primary operations on the RTE include ray intersection and 
ray reflection and transmission. The physical structure of the RTE is modeled as a plane with an 
aperture. Usually the aperture is circular, so the RTE looks like a circle of zero thickness in space. 
In local coordinates the RTE is represented as a plane equation with a separate test for the aperture. 
The equation of the planar object In local coordinates is simply 

2=0 

(1) 

The intersection of the RTE with a ray is straightforward. Let a ray in local coordinates be defined 
as In (2). 
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Then the intersection of the ray with the plane may be found, as in equation (3). 
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where z d = 0. Given the location of the planar Intersection (x,y, 0), a final test is made to ensure the 
point is inside the aperture region. 
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RTE Construction 

In the preferred embodiment, the construction of the RTE is described in local coordinates. In 
practice, rays\defined in world coordinates are first transformed into local coordinates using 
homogeneous transformation matrices, collectively called a world-to local transformation matrix. 
Likewise, resulting rays in local coordinates are transformed to world coordinates using an inverse 
of the previous homogeneous transformation matrix, called local-to-world transformation matrix, 

as would be well known to those skilled in the art ( G o o f o i e xample, 01jgs^A.S ., ^Intro trrjcT^n 

\ ^ 
io- Ray Tracing, AoadorAio Press; 1980 , for a discussion of homogeneous transformation matrices.) 

Since the goal of the RTE is to transform an input ray into an output ray for a given intersection 

point, we simply construct a mapping for a set of input/output ray pairs. Suppose an incident ray I 

and a transmitted ray rdefineain local Coordinates is given. The intersection point (x,y, 0) is denoted 

S. The plane that contains the op\ical axis (z-axis) and the point S is called the meridional plane. The 

plane parallel to the z-axis and perpendicular to the meridion at plane is referred to as the transverse 

plane. These planes are illustrated in Figure 3. The strategy for finding the mapping of I to T is to 

find the change in elevation and azimuth angles for the two rays with respect to the meridional and 

transverse pianos, so that I can be rotated into T. The first step is to compute a meridional plan 

coordinate basis matrix. For this coordinate system, X points from S to the origin. Y is 90 deg 

counter clockwise to X, and Z' is parallel tome original z-axis. These basis vectors are illustrated 

in Figure 4. \ 

For simplicity of presentation, in the following the prime notation is dropped for the 

meridional basis vectors X, Y, and Z'. Given S weVan compute the meridional coordinate basis 



matrix by finding the X, Y, and Z vectors as in equation (4). 
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If the norm of S is zero, the intersection point is at the origin, i.e., on the optical axis. In this case X 
in equation (4) is taken as being along the x-axis. Using this the meridional coordinate basis matrix 
M can be written as In (5). 



M = 



Xx y, z. 

X2 Yi Z2 
X } 7b Zb 



(5) 



To compute the desired rotation angles between / and T, first rotate the vectors using the meridional 
coordinate basis Mas in (6). 



T M = M r T 



(6) 



In (6) and in the following discussion, the vectors / and / are first normalized to unit length. Now, 
can be computed the rotation angle for each of these vectors to the meridional coordinate x-axis. 
Also can be computed the rotation angle between the XZ plane and each of these vectors. These 
calculations are made using equations (7) and (8), respectively. 
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0a = tan-i 
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4 A, 
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(7) 



a 



a = tan- 1 



x + 4 ; 



(8) 



In equations (7) and (8), A stands for either the vector / or T. The arc tangent function 
indicated in these equations is conveniently computed using the a atan2 function in C++ to 
automatically handle all possible values of numerator and denominator around the unit circle 
including the case of denominator equal to zero. After these two elevation angles and two azimuth 
angles have been computed, the difference can be computer as indicated in (9). These are the angles 
stored to map an incident ray direction to the transferred ray direction. 



d9 — 6 j — 6 T 

(9) 

da = a, - a T 

To summarize, for a given intersection point S, compute and save the elevation and azimuth 
angles d9 and da and do required to rotate the incident vector I into the direction of transferred vector 
T. All of the forgoing calculations would be known to those skilled in the art and for that reason are 
not described in detail. Discussed next, is the use of the RTE to map an input vector to an output 
vector. This illustrates the case of forward ray tracing. The case of reverse ray tracing is similarly 
computed, as would be known to one skilled in the art and is omitted for brevity. 
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RTE Optical Ray Tracing 

To perform an optical ray tracing of the RTE, first find where an incident vector I is 
intersected with the RTE at location S. This intersection point is computed using equation (3). Given 
this intersection point, the transmitted vector T is easily computed. In general, the intersection point 
S will not correspond to a previous constraint point, so some type of interpolation will be required. 
Two strategies to accomplish this are (1 ) to compute a least squares approximation function fit (e.g., 
2D polynomial or 2D B-Spline) for each of the azimuth and elevation angles or (2) to perform a 
Delaunay ^angulation of the intersection points in the XY plane and interpolate the points to obtain 
the angles. As would be known to those skilled in the art..', Suitable software may be developed for 
performing these calculations on a general purpose computer. The computation of 7 1 (which produces 
a ray tracing of the RTE) is performed as follows: 

1 . Intersect I with the RTE to find S 

2. Apply the transformations indicated in (4) to (6) and then rotate / by the azimuth and 
elevation angles associated with S. This rotated vector represents T in meridional coordinates. 

3. Apply the inverse coordinate transformation to T. 

The adjustment angles d6 and da are applied to the angles interpolated in step 2 of the 
algorithm. This provides a method to adjust the optical model based on actual outcomes. In a similar 
way, a paraxial power value could be assigned to the projection of the / end T rays onto the 
meridional and transverse planes. This alternate scheme has the benefit of describing the bending 
of rays in units of diopters that are familiar to visual optics professionals. 

12 



Strategy for RTE Updating Based on Historical Cases 

A strategy for updating the RTE portion of the optical model is based on the evaluation of 
historical cases. In a preferred embodiment, as described in this application, first examine the 
historical database for similar cases to a given present case under consideration. This similarity is 
measured in terms of both categorical data and continuous data.. Categorical data such as surgeon, 
sex, intervention procedure, etc., should be matched closely by all historical cases to be similar to 
the current case. This is accomplished by filtering the database for all records that match the same 
categories. Continuous data such as preoperative spectacle correction, corneal thickness, etc., is 
measured using a weighted Euclidean distance. Each of the 5 continuous data is normalized so that 
the mean and standard deviation of each random variable is zero and one, respectively . Only the N 
closest cases from the historical database are used to compute the RTE adaptation parameters. A 
reasonable value for the integer N may be determined during a first evaluation (described below). 

The RTE update algorithm is then: 

1 . Establish the categorical data for the current case under consideration and identify 
the continuous data. 

2. Filter the historical database so that we only consider cases in the same category. 

3. For all filtered records, find the N closest records using normalized distances. 

4. For these N closest records, find the average back-calculated RTE for the actual 
postoperative wavefront aberration (See RTE Arithmetic Operations below.) 

5. Use the average of the RTE deltas in the Optical Modeling of the current case 
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I It is often possible to have the case of time-varying parameters, that is, cases where the 
statistics slowly vary over time. To account for this in our calculations of the updated RTE, a 
weighting may be applied to the average calculation ab that older cases are not weighted as much 
as more recent cases. 

RTE Arithmetic Operations 

RTE arithmetic operations that help facilitate AVIA calculations are: 

Addition, subtraction, multiplication, and division of two RTEs. 

Addition, subtraction, multiplication, and division of an RTE and a scalar. 

These operations (addition, subtraction, multiplication, and division) on two RTEs are 
performed in a point-wise fashion. In general, the stored sample points from two RTEs will not be 
in the same locations, so one or both of the sample point arrays must be interpolated to obtain points 
at matching locations. 

These operations (addition, subtraction, multiplication, and division) on a RTE and a scalar 
are also performed in a point- wise fashion. In this case there is no need to interpolate to align sample 
points as the same scalar operation is applied to all points in the RTE. 

Implementation of the RTE using the C++ programming language permits these operations 
to be specified directly using standard symbols of "+", "- . and "7". This makes the generation and 
maintenance of RTE aware programs very efficient. For example, the mean of four RTEs can be 
expressed as: 

(R1+ R2+ R3+ R4) 
MR = (10) 
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1 Using an appropriate class implementation of a RTE, this expression is a valid C++ statement 

2 and its intent is easily understood. 
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